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Abstract 

The hadron mass spectrum is calculated in lattice QCD using a novel fat- 
link clover fermion action in which only the irrelevant operators in the fermion 
action are constructed using smeared links. The simulations are performed on 
a 16^ X 32 lattice with a lattice spacing of a = 0.125 fm. We compare actions 
with n = 4 and 12 smearing sweeps with a smearing fraction of 0.7. The n = 4 
Fat-Link Irrelevant Clover (FLIC) action provides scaling which is superior 
to mean-field improvement, and offers advantages over nonperturbative 0{a) 
improvement, including a reduced exceptional configuration problem. 
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I. INTRODUCTION 



The origin of the masses of hght hadrons represents one of the most fundamental chal- 
lenges to the theory of strong interactions, Quantum Chromodynamics (QCD). Despite the 
universal acceptance of QCD as the basis from which to derive hadronic properties, there 
has been slow progress in understanding the generation of hadron mass from first principles. 
Solving the problem of the hadronic mass spectrum would allow considerable improvement 
in our understanding of the nonperturbative nature of QCD. The only available method at 
present to derive hadron masses directly from QCD is a numerical calculation on the lattice. 
In the last few years impressive progress has been made both in computer hardware and in 
developing more efficient algorithms, bringing realistic simulations of hadronic observables 
with sufficiently large volumes, small quark masses and fine enough lattices within reach. 

The high computational cost required to perform accurate lattice calculations at small 
lattice spacings, however, has led to increased interest in quark action improvement. In order 
to avoid the famous doubling problem, Wilson originally introduced an irrelevant (energy) 
dimension-five operator (the "Wilson term" ) to the standard "naive" lattice fermion action, 
which explicitly breaks chiral symmetry at 0{a). To extrapolate reliably to the continuum, 
simulations must be performed on fine lattices, which are therefore very computationally 
expensive. The scaling properties of the Wilson action at finite a can be improved by 
introducing any number of irrelevant operators of increasing dimension which vanish in the 
continuum limit. 

The Sheikholeslami-Wohlert (clover) action introduces an additional irrelevant 
dimension-five operator to the standard Wilson quark action, 

Ssw = Sw- ^^'^^^^ ilj{x)a^^Ff,^^lj{x) , (1) 
where S'w is the standard Wilson action. 



5'w = ip{x) 



5Z ( 7m - ^raA^ ) + m 



^(x) , (2) 



and are the standard covariant first and second order lattice derivatives, 

= ^ [Uf,{x)'ilj{x + /i) - Ul{x - fi)ij{x - fi) 
A^ip{x) = \ Uf,{x)ip{x + /i) + Ul{x - ij)ipix - /i) - 2ip{x) 

(X 

and Csw is the clover coefficient which can be tuned to remove 0{a) artifacts, 

^ _ 1 1 at tree- level , 
sw I mean-field improved , 

with Uo the tadpole improvement factor which corrects for the quantum renormalization 
of the operators. Nonperturbative (NP) 0{a) improvement tunes Csw to all powers in 
g"^ and displays excellent scaling, as shown by Edwards et al. [0], who studied the scaling 
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properties of the nucleon and vector meson masses for various lattice spacings (see also 
Section |^ below). In particular, the linear behavior of the NP-improved clover actions, 
when plotted against a?, demonstrates that 0{a) errors are removed. It was also found in 
Ref. [§1 that a linear extrapolation of the mean-field improved data fails, indicating that 
0{a) errors are still present. 

A drawback to the clover action, however, is the associated problem of exceptional con- 
figurations, where the quark propagator encounters singular behavior as the quark mass 
becomes small. In practice, this prevents the use of coarse lattices (/3 > 5.7 ~ a ^ 0.18 fm) 
Furthermore, the plaquette version of F^^, which is commonly used in Eq. (|Tp, has 
large 0{a^) errors, which can lead to errors of the order of 10% in the topological charge 
even on very smooth configurations 0. 

The idea of using fat links in fermion actions was first explored by the MIT group ^ and 
more recently has been studied by DeGrand et al. ||,|,|l3], who showed that the exceptional 



configuration problem can be overcome by using a fat-link (FL) clover action. Moreover, the 
renormalization of the coefficients of action improvement terms is small. In principle it is 
acceptable to smear the links of the relevant operators. The symmetry of the APE smearing 
process ensures that effects are 0{a^). The factors multiplying the link and staple ensure the 
leading order term is 6*"^"^^ , an element of SU(3). Issues of projecting the smeared links to 
SU(3) are 0{a?) effects and therefore correspond to irrelevant operators [11|. However, the 
net effect of APE smearing the links of the relevant operators is to remove gluon interactions 
at the scale of the cutoff. While this has some tremendous benefits, the short- distance quark 
interactions are lost. As a result decay constants, which are sensitive to the wave function 
at the origin, are suppressed. 

A possible solution to this is to work with two sets of links in the fermion action. In the 
relevant dimension-four operators, one works with the untouched links generated via Monte 
Carlo methods, while the smeared fat links are introduced only in the higher dimension 
irrelevant operators. The effect this has on decay constants is under investigation and will 
be reported elsewhere. 

In this paper we present the first results of simulations of the spectrum of light mesons 
and baryons using this variation on the clover action. In particular, we will start with the 
standard clover action and replace the links in the irrelevant operators with APE smeared 
T^ , or fat links. We shall refer to this action as the Fat-Link Irrelevant Clover (FLIC) 



action. Although the idea of using fat links only in the irrelevant operators of the fermion 
action was developed here independently, suggestions have appeared previously . To the 



best of our knowledge, this is the first report of lattice QCD calculations using this novel 
fermion action. 

In Section H we describe the gauge action used in our lattice simulations, while Section |T| 
contains the procedure for creating the FLIC fermion action. Our results are presented in 
Section 0, and finally in Section |V] we draw some conclusions and discuss future work. 

II. THE GAUGE ACTION 

The simulations are performed using a mean-field improved, plaquette plus rectangle, 
gauge action on a 16^ x 32 lattice at (3 = = 4.60, providing a lattice spacing a = 
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0.125(2) fm determined from the string tension with ^Ja = 440 MeV. The tree-level 0(a^)- 
Symanzik-improved gauge action [14| is defined as 



= ^E^e tr(l - U^^ix)) -J^Y.'^e tr(l - U,,,,{x)) , (4) 

sq ^^^0 rcct 

where the operators f/sq(x) and Urectix) are defined as 

Usqix) = U^{x)U,{x + fi)Ul{x + i))Ul{x) , (5a) 
f^rcct(a;) = U^,{x)U^{x + fi)U^{x + U + fl) 
xUl{x + 2i))Ul{x + i))Ul{x) 
+ U^{x)U^{x + fi)Uy{x + 2/i) 

xUl{x + fi + v)Ul{x + v)Ul{x) . (5b) 

The link product f/rect(a;) denotes the rectangular 1x2 and 2x1 plaquettes, and for the 
tadpole improvement factor we employ the plaquette measure 

/I \^'^ 
Mo= (^-7^etr(f/,q)j . (6) 

Gauge configurations are generated using the Cabibbo-Marinari pseudoheat-bath algorithm 
with three diagonal SU(2) subgroups looped over twice. Simulations are performed using a 
parallel algorithm with appropriate link partitioning |T^. A total of 50 configurations are 



used in this analysis, and the error analysis is performed by a third-order, single-elimination 
jackknife, with the P^r degree of freedom (A^df) obtained via covariance matrix fits. 

III. FAT-LINK IRRELEVANT FERMION ACTION 

Fat links HH] are created by averaging or smearing links on the lattice with their nearest 



neighbours in a gauge covariant manner (APE smearing). The smearing procedure |jT2 
replaces a link, U^{x), with a sum of the link and a times its staples 
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U^.{x) ^ f/;(x) = {l-a)U^{x) + '^Y.[UAx)U^.{x + ^o)Ul{x + ^ia) 

^ i/=i 

+ Ul{x — va)U^{x — ua)Uy{x — ua + jJLa) , (7) 
followed by projection back to SU(3). We select the unitary matrix Uj^^ which maximizes 

^etr(f/JLf/;), 

by iterating over the three diagonal SU(2) subgroups of SU(3). We repeat the combined 
procedure of smearing and projection n times. We create our fat links by setting a = 0.7 
and comparing n = 4 and 12 smearing sweeps. The mean-field improved FLIC action now 
becomes 
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C-FL 



qFL 



(8) 



where F^j, is constructed using fat links, Mq is calculated in an analogous way to Eq 
and where the mean-field improved Fat-Link Irrelevant Wilson action is 



C-FL 



X X,fl 

jjFL 



7m 



Uo 



'4>{x + fi) 



— r 



,,FL 



,,FL 



ij{x - fi) 



ip{x — fl) 



(9) 



with K = l/(2m -|- 8r). We take the standard value r = 1. Our notation uses the Pauli 
representation of the Dirac 7-matrices defined in Appendix B of Sakurai |T^. In particular, 
the 7-matrices are hermitian and a^^ = [7^, 7i,]/(2i). 

As reported in Table |, the mean-field improvement parameter for the fat links is very 
close to 1. Hence, the mean-field improved coefficient for Csw is expected to be adequate[|. In 
addition, actions with many irrelevant operators {e.g. the D234 action) can now be handled 
with confidence as tree-level knowledge of the improvement coefficients should be sufficient. 
Another advantage is that one can now use highly improved definitions of (involving 



terms up to Uq), which give impressive near-integer results for the topological charge JT 

In particular, we employ an O(a^) improved definition of F^jy in which the standard 
clover-sum of four 1x1 Wilson loops lying in the /i, u plane is combined with 2x2 and 



3x3 Wilson loop clovers. Bilson-Thompson et al. [0] find 



p = 

8 




20m^ 



W 



2x2 



+ 



W 



3x3 



h.c. 



(10) 



Traceless 



where W"'^^ is the clover-sum of four n x n Wilson loops and F^j^^, is made traceless by 
subtracting 1/3 of the trace from each diagonal element of the 3x3 color matrix. This 
definition reproduces the continuum limit with 0{a^) errors. On approximately self-dual 
configurations, this operator produces integer topological charge to better than 4 parts in 
10^ 

Work by DeForcrand et al. []TB[ suggests that 7 cooling sweeps are required to approach 
topological charge within 1% of integer value. This is approximately 16 APE smearing 
sweeps at a = 0.7 ||T^. However, achieving integer topological charge is not necessary for 
the purposes of studying hadron masses, as has been well established. To reach integer topo- 
logical charge, even with improved definitions of the topological charge operator, requires 
significant smoothing and associated loss of short- distance information. Instead, we regard 
this as an upper limit on the number of smearing sweeps. 

Using unimproved gauge fields and an unimproved topological charge operator. Bonnet et 
al. found that the topological charge settles down after about 10 sweeps of APE smearing 



Our experience with topological charge operators suggests that it is advantageous to include uq 
factors, even as they approach 1. 
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n 









0.88894473 


0.62445197 


4 


0.99658530 


0.98641100 


12 


0.99927343 


0.99709689 



TABLE I. The value of the mean link for different numbers of APE smearing sweeps, n, at 
a = 0.7. 



at a = 0.7. Consequently, we create fat links with APE smearing parameters n = 12 
and a = 0.7. This corresponds to ~ 2.5 times the smearing used in Refs. |]^,|^. Further 
investigation reveals that improved gauge fields with a small lattice spacing (a = 0.125 fm) 
are smooth after only 4 sweeps. Hence, we perform calculations with 4 sweeps of smearing 
at a = 0.7 and consider n = 12 as a second reference. Table | lists the values of Uq^ for 
n = 0, 4 and 12 smearing sweeps. 

We also compare our results with the standard Mean-Field Improved Clover (MFIC) 
action. We mean- field improve as defined in Eqs. ^ and ^ but with thin links throughout. 
The standard Wilson-loop definition of F^^ is used. 

A fixed boundary condition is used for the fermions by setting 

Ut{x,nt) = and Uf^{x,nt)=0 Vx (11) 

in the hopping terms of the fermion action. The fermion source is centered at the space-time 
location {x, y, z, t) = (1, 1, 1, 3), which allows for two steps backward in time without loss of 
signal. Gauge-invariant gaussian smearing in the spatial dimensions is applied at the source 
to increase the overlap of the interpolating operators with the ground states. 

IV. RESULTS 

Hadron masses are extracted from the Euclidean time dependence of the calculated two- 
point correlation functions. For baryons, the correlation functions are given by 

(G(t;p,r)) = J2e-'^^-^ r^"(fi|T[x"(x)x^(0)]|fi) , (12) 

X 

where x are standard baryon interpolating fields, Q represents the QCD vacuum, F is a 4 x 4 
matrix in Dirac space, and a, (3 are Dirac indices. At large Euclidean times one has 

(G(t; pj))^^ e-^''* tr [F(-Z7 ■ p + M)] , (13) 

where Z represents the coupling strength of x(0) to the baryon, and Ep = (p^ + M^)^/^ is 
the energy. Selecting p = and F = (1 + 74) /4, the effective baryon mass is then given by 
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FIG. 1. Effective mass plot for the nucleoli for the FLIC action with 4 sweeps of smearing 
at Q = 0.7 from 200 configurations. The five sets of points correspond to the k values listed in 
Table 0, with K increasing from top down. 



M{t + 1/2) = log[G(t)] - \og[G{t + 1)] . (14) 

Meson masses are determined via analogous standard procedures. The critical value of k, Kc, 
is determined by linearly extrapolating as a function of rriq to zero. The quark masses 
are defined by = (I/k — 1/kc) /(2a), and the strange quark mass was taken to be the 
second heaviest quark mass in each case. 

Figure ^ shows the nucleon effective mass plot for the FLIC action when 4 APE smearing 
sweeps at a = 0.7 are performed on the fat links ("FLIC4"). The effective mass plots for 
the other hadrons are similar, and all display acceptable plateau behavior. Good values of 
X^/Ndp are obtained for many different time-fitting intervals as long as one fits after time 
slice 8. All fits for this action are therefore performed on time slices 9 through 14. For the 
Wilson action and the FLIC action with n = 12 ("FLIC12") the fitting regimes used are 
9-13 and 9-14, respectively. 



The values of k used in the simulations for all quark actions are give in Table |T|. We 
have also provided the values of Kcr for these fermion actions when using our mean-field 
improved, plaquette plus rectangle, gauge action at /5 = 4.60. We have mean- field improved 
our fermion actions so we expect the values for k^t to be close to the tree-level value of 0.125. 
Improved chiral properties are seen for the FLIC and MFIC actions, with FLIC4 performing 
better than FLIC12. 

The behavior of the p, nucleon and A masses as a function of squared pion mass is 
shown in Fig. |^ for the various actions. The first feature to note is the excellent agreement 
between the FLIC4 and FLIC12 actions. On the other hand, the Wilson action appears 
to lie somewhat low in comparison. It is also reassuring that all actions give the correct 



7 



K ttItt a nip a ttt-n a ttt-a a 

0.1260 0.5797(23) 0.7278(39) 1.0995(58) 1.1869(104) 

0.1266 0.5331(24) 0.6951(45) 1.0419(64) 1.1387(121) 

0.1273 0.4744(27) 0.6565(54) 0.9709(72) 1.0816(152) 

0.1279 0.4185(30) 0.6229(65) 0.9055(82) 1.0310(194) 

0.1286 0.3429(37) 0.5843(97) 0.8220(102) 0.9703(286) 



TABLE II. Values of k and the corresponding vr, p, N and A masses for the FLIC actfon with 4 
sweeps of smearing at a = 0.7. The value for Hcr is provided in Table ID. A string tension analysis 
provides a = 0.125(2) fm for ^ = 440 MeV. 





Wilson 


FLIC12 


FLIC4 


MFIC 




0.1346 


0.1286 


0.1260 


0.1196 


K2 


0.1353 


0.1292 


0.1266 


0.1201 


K3 


0.1360 


0.1299 


0.1273 


0.1206 


K4 


0.1367 


0.1305 


0.1279 


0.1211 


K5 


0.1374 


0.1312 


0.1286 


0.1216 




0.1390 


0.1328 


0.1300 


0.1226 



TABLE III. Values of k and Krr for the four different actions. 



mass ordering in the spectrum. The value of the squared pion mass at 771^^/ nip = 0.7 is 
plotted on the abscissa for the three actions as a reference point. This point is chosen in 
order to allow comparison of different results by interpolating them to a common value of 
mTr/rrip = 0.7, rather than extrapolating them to smaller quark masses, which is subject to 
larger systematic and statistical uncertainties. 

The scaling behavior of the different actions is illustrated in Fig. ^. The present results 
for the Wilson action agree with those of Ref . 0] . The first feature to observe in Fig. ^ is 
that actions with fat-link irrelevant operators perform extremely well. For both the vector 
meson and the nucleon, the FLIC actions perform significantly better than the mean-field 
improved clover action. It is also clear that the FLIC4 action performs systematically better 
than the FLIC12. This suggests that 12 smearing sweeps removes too much short-distance 
information from the gauge-field configurations. On the other hand, 4 sweeps of smearing 
combined with our O(a^) improved Fp^ provides excellent results, without the fine tuning 
of Csw in the NP improvement program. 

Notice that for the p meson, a linear extrapolation of previous mean-field improved clover 
results in Fig. ^ passes through our mean-field improved clover result at a^a ~ 0.08 which 



8 



2.0 



> 



CO 



A 

LA 

N 


1 


1^ 


1 


* * 


1 


1 

* 

m 




P 


1 




IS 






• FLIC12 






n 


1 


* 
■ 


1 


■ FLIC4 
° Wilson 





0.0 0.2 0.4 0.6 0.8 1.0 

(GeV^) 

FIG. 2. Masses of the nucleon, A and p meson versus for the FLIC4, FLIC12 and Wilson 
actions. 



lies systematically low relative to the FLIC actions. However, a linear extrapolation does 
not pass through the continuum limit result, thus confirming the presence of significant 0{a) 
errors in the mean-field improved clover fermion action. While there are no NP-improved 
clover plus improved glue simulation results at a'^a ~ 0.08, the simulation results that are 
available indicate that the fat-link results also compete well with those obtained with a 
NP-improved clover fermion action. 

Having determined FLIC4 is the preferred action, we have increased the number of 
configurations to 200 for this action. As expected, the error bars are halved and the central 
values for the FLIC4 points move to the upper end of the error bars on the 50 configuration 
result, further supporting the promise of excellent scaling. 

Finally, in Fig. ^ we compare the convergence rates of the different actions by plotting the 
number of stabilized biconjugate gradient iterations required to invert the fermion matrix 
as a function of mj^/ rrip. For any particular value of rrij^/ nip, the FLIC actions converge faster 
than both the Wilson and mean-field improved clover fermion actions. Also, the slopes of 
the FLIC lines are smaller in magnitude than those for Wilson and mean-field improved 
clover actions, which provides great promise for performing cost effective simulations at 
quark masses closer to the physical values. Problems with exceptional configurations have 
prevented such simulations in the past. 
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cr 

FIG. 3. Nucleoli and vector meson masses for the Wilson, NP-improved, mean-field clover and 
FLIC actions. Results from the present simulations, indicated by the solid points, are obtained by 
interpolating the results of Fig. ^ to mT^/nip = 0.7. The fat links are constructed with n = 4 (solid 
squares) and n = 12 (stars) smearing sweeps at a = 0.7. The solid triangles are results for the 
FLIC4 action when 200 configurations are used in the analysis. The FLIC results are offset from 
the central value for clarity. Our MF clover result at a^a ~ 0.08 lies systematically low relative to 
the FLIC actions. 

V. CONCLUSIONS 

We have examined the hadron mass spectrum using a novel Fat-Link Irrelevant Clover 
(FLIC) fermion action, in which only the irrelevant, higher-dimension operators involve 
smeared links. One of the main conclusions of this work is that the use of fat links in the 
irrelevant operators provides excellent results. Fat links promise improved scaling behavior 
over mean- field improvement. This technique also solves a significant problem with 0{a) 
nonperturbative improvement on mean field-improved gluon configurations. Simulations 
are possible and the results are competitive with nonperturbative-improved clover results 
on plaquette-action gluon configurations. We have found that minimal smearing holds the 
promise of better scaling behavior. Our results suggest that too much smearing removes 
relevant information from the gauge fields, leading to a poorer performance. Fermion ma- 
trix inversion for FLIC actions is more efficient and results show no sign of exceptional 
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FIG. 4. Average number of stabilized biconjugate gradient iterations for the Wilson, FLIC and 
mean-field improved clover (MFIC) actions plotted against m-,j/mp. The fat links are constructed 
with n = A (solid squares) and n = 12 (stars) smearing sweeps at q = 0.7. 



configuration problems. 

This work paves the way for promising future studies. It will be of great interest to 
consider different lattice spacings to further test the scaling of the fat-link actions. Further- 
more, the exceptional configuration issue can be explored by pushing the quark mass down 
to lower values. We are currently performing simulations at m^/mp = 0.36 and these results 
will available soon. A study of the spectrum of excited hadrons using the FLIC actions is 
also currently in progress pi] . 
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